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GBM Detectors 


Nal Detectors. The Nal(TI) Scintalators are 127 
mm in diameter and 12.7 mm thick. The 
detectors covers the 8 keV - 1 MeV band. 



Pulsed Frequency and Pulsed Flux Monitoring 


The analysis of GBM observations of pulsars presents two main 
challenges: 

•The background rates are normally much larger than the source 
rates, and have large variations. 

•The responses of the detectors to a source are continuously 
changing because of Fermi’ s continuous re-orientation. 

The initial steps of the pulsed data analysis are: 

•Data Screening 

•Background subtraction of the Nal detector count rates 
•Determination of fluxes from remaining rates 


Mark H. Finger 


MAXI Workshop 


1 December 2010 
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Background Subtraction 



The rates in each 
channel of the 12 Nal 
detectors are fit with a 
model with the following 
components: 

•Models for bright 
sources. 

•A stiff empirical model 
that contains the low- 
frequency component of 
the remaining rates. 

The fits are made 
independently for each 
channel and subtracted 
from the rates. 


Mark H. Finger 


MAXI Workshop 


1 December 2010 



Estimating Pulsed Fluxes 


For a given source we combine the rate residuals over detectors 
and obtain an estimate of the variable part of the source flux. 
Using a model of the source spectrum and the time dependent 
detector responses we compute the source induced rate ^ 
expected in detector / at time tk if the source has unit flux in the 
channel’ s energy range. The variable part of the flux/fcis then 
estimate by minimizing 


where fy s the residual rates and cr|| le associated errors. 




Pulse Searches 


We have implemented two different pulse search strategies: 

•Daily Blind Search. For this we compute fluxes from a days data 
for 24 source directions equally spaced on the galactic plane. For 
each direction we do an FFT based search from 1 mFIz to 2 Hz. 

•Source Specific Searches. These are searches over small ranges 
of frequency and sometimes frequency rate based on phase 
shifting and summing pulse profiles that are made from short 
intervals of data, using barycentered and possibly orbitally 
corrected times. Typical exposure times are ~40 ks/day. 



Corrected Power 


Blind Pulse Search 



Blind pulse search in 20-50 keV band, for 2010 January 8. 





http://gammaray.nsstc.nasa.gov/gbm/science/pulsars 
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GBM Pulsar Project 


GBM-related pulsar publications from the Gamma-Ray Astronomy team are listed here . 


GBM Accreting Pulsar Histories 

For each source wc plot the history of pulse frequency and pulsed flux measured using the Fermi Gamma-Ray Burst Monitor (GBM) Nal detectors. For these measurements we use 
the CTIME data which normally has 0.256 s time bins, and eight energy channels. Our analysis normally uses channels 1 (12-25 kcV) and 2 (25-50 keV). The integration interval 
| used varies from source to source, ranging from one to four days. For eclipsing systems each egress to ingress interval is divided into an integral number of equal parts, with no 
measurement made during the eclipse. The measured frequencies arc b ary centered. For sources where the binary orbit is known, the frequencies arc corrected for the binary motion. 

I The R.M.S. pulsed flux is given in the energy band in which the pulse search was made. This usually includes only the first and second harmonics. These results arc preliminary. 
Please contact Mark Finger for further information. 

Detected Persistent Pulsars 


ISource Name 
|& GBM link 

LB (deg) 

BII (deg) 

Period (s) 

Swift link 

RXTE link 

MAXI link 

Simbad link 

|GX 1+4 

1.94 

4.79| 

159.7 

GX 1-14 

gxl+4 

[GX 1+4 

V- V21 16 Oph 

[Her X-l 

58.20 

37.50 

1.24 

Her X-l 

hcrxl 

|HcrX-l 

V* HZ Her 

jvcla X-l 

263.06 

3.90 

283.5 

Vela X-l 

vclaxl 

[Vela X-l 

V* GP Vcl 

|Ccn X-3 

292.10; 

0.30 

4.80 

Cen X-3 

ccnx3 

|Ccn X-3 

V* V779 Cen 

|GX 301-2 

300. IB 

1.25! 

681.6 

GX 301-2 

gx301-2 

[GX 301-2 

V*BPCru 

|4U 1626-67 

321.79 

-13.09 

7.67 

4U 1626-67 

Kl«7-$73 

|4U 1626-67 

V* KZ TrA 

|4U 1538-52 

327.42 

2.16| 

525.0 

H 1538-522 

X1538-522 

|4U 1538-52 

V* OV Nor 

|QAQ 1657-415 

344.40 

0.3 l| 

37.1 

EXO 1657-419; 

X1657-415 

joAQ 1657-415 

OAO 1657-41 


Detected Transient Pulsars 


Source Name 
& GBM Link 

LB (deg) 

BB (deg) 

Period (s) 

Swift link 

RXTE link 

MAXI link 

Simbad link 

IGR J18179-1621 

14.60 

-0.22 

11.82 

IGR J18179-1621 




XTE J 1858+034 

36.82 

-0.06 

219 

XTE J 1858+034 

xtei 1858+034 

XTE J 1858+034 

fKRL2007bl 348 

4U 1901+03 

37.16 

-1.25 

2.76 

4U 1901+03 

X1901+031 

4U 1901+03 

4U 1901+03 

IGR J 1 9294+ 1816 

53.44 

0.20 

12.45 

IGR J19294+1816 

igri 1 9294+ 1816 


IGR J 1 9294+ 1816 

XTE J 1946+274 

63.21 

1.40 

15.8 

XTE J 1946+274 

x 1942+274 

XTE J 1946+274 

3 A 1942+274 

EXO 2030+375 

77.15 

-1.24 

41.33 

EXO 2030+375 

cxo2030+375 

EXO 2030+375 

V* V2246 Cvg 

SAX J2 103 .5+4545 

87.13 

-0.68 

358.61 

SAX J2 103 .5+4545 

sax i2 103 .5+4545 

SAX J2 103 .5+4545 

GSC 03588-00834 

CcpX-4 

99.01 

3.31 

66.32 

Ginca 2138+56 

cepx4 

CcdX-4 

V* V490 Ccd 

4U 01 15+63 

125.92 

1.03 

3.61 

4U 01 15+634 

xO 115+634 

4U 01 15+63 

V* V635 Cas 

V 0332+53 

146.05 

-2.19 

4.37 

V 0332+53 

V0332+53 

V 0332+53 

V* BO Cam 

RX JQ440 .9+4431 

159.85 

-1.27 

205.26 

LSV+44 17 

lsv+4417 

LSV+4417 

TYC 2905-121-1 

A 0535+26 

181.50 

-2.64 

103.5 

1 A 0535+262 

X0535+262 

A 0535+262 

HIP 26566 

MXB 0656-072 

220.20 

-1.76 

160.4 

MXB 0656-072 

X0656-072 

MXB 0656-072 

2E 0655.8-0708 

GS 0834-430 

262.02 

-1.51 

12.3 



GS 0834-430 






Spin Frequency (mHz) 



Pspin = 283 s, Porbit = 8.96 d, wind-fed supergiant 



GBM Earth Occultation Monitoring 

The observed count rate in a 
240-s window of data is fitted 
by a model consisting of a 
quadratic background plus 
source terms. 


r (t , E C h) - bo(E C h) + bi(E C h) (t - to) + 

n 

fe(Ech) (t- t 0 ) 2 + J>(E ch ) Si(t, Ech) (2) 


The source count rate is modeled using an assumed energy spectrum, 
convolved with the atmospheric transmission function and folded 
through the detector response matrix 


S(t,E ch )= R(E ph ,E ch ,t) |T(Ep h ,t)) [ f (E ph )dE ph 

V -'Eph 

(3) 








GBM Earth occultation Monitoring (2) 


• Using the first 3 years of GBM data, a catalog of 
209 sources was monitored, with 99 sources 
detected, including 41 of 52 LMXBs, 31 of 39 
HMXBs, 12 of 19 BHCs, 12 of 71 AGN, the 
Ophiuchus Cluster, the Sun, and the Crab. 

• 7 BHCs, 1 AGN (Cen A), and the Crab are 
detected above 100 keV. 

• Typical source exposure times are 3ks/day. 


Wilson-Hodge et al. 2012, ApJ, Accepted, arXiv:1201.3585 



http://heastro.phys.lsu.edu/gbm 









Persistent Sources Detected 


Name 

Spin Period (s) Orbital Period (d) 

Type 

Her X-1 

1.24 

1.70 

Disk-fed LMXB (Eclipsing) 

Cen X-3 

4.80 

2.09 

Disk-fed Supergiant (Eclipsing) 

4U 1626-67 

7.63 

0.023 

Disk-fed LMXB (Super- 
Compact) 

OAO 1657-415 

37.1 

10.4 

Wind-fed Supergiant (Eclipsing) 

GX 1+4 

158 

1161 

Symbiotic XRB (red giant+ns) 

Vela X-1 

283 

8.96 

Wind-fed Supergiant (Eclipsing) 

4U 1538-52 

525 

3.73 

Wind-fed Supergiant (Eclipsing) 

GX 301-2 

686 

41.5 

Wind-fed Supergiant 






Cen X-3 



Time (MJD) 


spin = 4.80 s, Porbit = 2.89 d, disk-fed supergiant 



1 2*50 keV Pulsed Flux 

(keV cm ? s ') Frequency (mHz) 


6.55 


GX 1+4 



spin - 158 s, Porbit _ 1161 d, Sy m biotic XRB 



Transient Sources Detected (1) 


Name 

Spin Period (s) 

Orbital Period (d) 

Type 

4U 1901+03 

2.76 

22.5 

Be/X-ray Binary 

4U 0115+634 

3.61 

24.3 

Be/X-ray Binary 

V 0332+53 

4.37 

34.2 

Be/X-ray Binary 

IGR J 1 81 79-1621 

11.82 



GS 0834-430 

12.3 


Be/X-ray Binary 

IGR J19294+1816 

* 

12.4 

117.2 

SFXT or Be/X-ray Binary 

XTEJ 1946+274 

15.8 

169.2 

Be/X-ray Binary 

2S 1417-624 

17.5 

42.1 

Be/X-ray Binary 

Swift J051 3.4-6547 

* 

27.3 

? 1 

likely Be/X-ray Binary (in 
LMC) 

EXO 2030+375 

41.3 

46.0 

Be/X-ray Binary 

Cep X-4 

66.3 

? 

Be/X-ray Binary 

GRO J 1008-57 

93.7 

248 

Be/X-ray Binary 

A 0535+26 

103 

111.1 

Be/X-ray Binary 




Transient Sources Detected (2) 


Name 

Spin Period (s) 

Orbital Period (d) 

Type 

MXB 0656-072 

160 

? 

Be/X-ray Binary 

LSV+44 1 7 

205 

-150 

Be/X-ray Binary 

XTEJ 1858+034 

219 


Be/X-ray Binary 

GX 304-1 

276 

132.5 

Be/X-ray Binary 

SAX 

J21 03.5+4545 

352 

12.68 

Be/X-ray Binary 

A 1118-615 

407 

? 

Be/X-ray Binary 

MAXI J 1409-6 19 

504 






Transient Outbursts with Detected Pulsations 





Transient Sources Detected with Earth Occultation 

m XTE J 1 946+274 
Z XTE J1 858 + 034 
W XTE J 1 752-223 
XTE J1701 -462 
• X Per 

• V0332+53 

Ser X— 1 
SWIFT J 1843.5-0343 
SWIFT J 1 842.5- 1 124 
SWIFT J1539. 2-6227 

m smc x-i 

A SAX 03.5 + 4545 

Mrk 421 
• MXB 0656-072 
MAXI J 1659- 152 
#MAXI J 1409-61 9 
LMC X — 4 
IGR J 1 8539 + 0727 
IGR J 1 7473-272 1 
IGR J 1 7464-32 1 3 
IGR J17091 -3624 
IGR J 1 6479-45 1 4 
IGR J 1 63 1 8-4848 
• Her X-1 

HT 1900.1-2455 
H 1 820-303 
HI 705-440 
HI 608-522 
H 1 4 1 7-624 
GX 354-0 
GX 339-4 

• GX 304-1 

• GX 301-2 

GX 1+4 
GS 1843+00 
GRS 1758-258 
GRO J 1 008-57 
GK Per 
GINGA 2138 + 56 
EXO 2030+375 
EXO 1745-248 
Cir X-1 
Cen X — 4 
Aql X-1 
AX J 1 631 . 9-4752 
A1 1 18-61 
A0535+262 
4U 1954+31 
4U 1901+03 
41) 1702-429 
4U 1636-536 
4U 1630-472 
4U 0614+09 
4U 01 15+634 
3A 01 14+650 
IE 1 148 1 -6141 
1 A 1742-294 

i i | 1 i i | i l I | 1 I 1 l m i -| i l 




■ Peak Flux 50-150 mCrab _ 

■ Peak Flux 150-500 mCrab — 

■ Peak Flux > 500 mCrab — 
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Torque Switching 






Flux / (HEAO-1 Flux) 


4U 1626-27 torque switching 


Year 

1980 1985 1990 1995 2000 2005 



Camero-Arranz et al. 2010 


2008 torque reversal 





4U 1626-27 torque switching 


Things to note: 

• The reversals are rapid compared to their separation. 

• During the spin-down the frequency rate increased while the 
flux decreased. This is inconsistent with a monotonic 
relationship between flux and frequency rate. 

• The spin-up to spin down reversal occurred at a higher flux 
than the spin-down to spin-up reversal. This is inconsistent 
with a single-valued relationship between flux and frequency 
rate. 



Cen X-3 torque switching 





Pulse Frequency (mHz) 


GX 1+4 Torque Switching 





GX 1+4 Torque Switching 


Things to note: 

• During the spin-down changes in frequency rate are anti- 
correlated with changes in flux. This supports the proposal a 
retrograde accretion disk during spin-down. 

• Two short-lived episodes of spin-up were seen with BATSE. 



OAO 1657-415 torque switching 


OAO 1657-415 
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GX 301-2 torque switching? 





Spin-down in wind accretion 

In spherically symmetric accretion the torque on the neutron 
star is estimated as (Lipunov 1992) 

r = 27 t/i> = —k/j, 2 r~g 

where k < 1 and the co-rotation radius is 

rco = 


Source 


OAO 1657-415 


B(10 12 G) period (s) 


37.7 


measured dv/dt calculated dv/dt 
(Hz/s) (Hz/s) 


-2.0 x 10‘ 12 


-5.9 x 10‘ 13 
x (B/10 13 ) 2 


GX 1+4 


158.4 


-4.4 x 10‘ 12 


-3.4 x10- 14 
x (B/10 13 ) 2 













Flux (10 8 erg cm' 7 s' 1 ) 


A0535+26 QPO 


2009 Dec 14 





Coherence Function y 2 Hard Lag (°) 


A0535+26 QPO 


2009 Dec 14 





Conclusions 

The full sky coverage of GBM enables long term monitoring of 
the brighter accreting pulsars allowing: 

•Precise measurements of spin frequencies and orbital 
parameters. 

•Study of spin-up or spin-down rates and hence the flow of 
angular momentum. 

•Detection and study of new transient sources or new outburst 
of known transients. 


GBM Pulsar Project 

http://gammaray.nsstc.nasa.gov/gbm/science/pulsars/ 



